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A  THEORETICAL  AND  EXPERIMENTAL  INVESTIGATION  OF  THE 
OXIDATION  OF  MOLYBDENUM  AT  TEMPERATURES 
AT  WHICH  ITS  TRIOXIDE  IS  VOLATILE 

By  David  R.  Schryer  and  Gerald  D.  Walberg 
Langley  Research  Center 

SUMMARY 


A  theoretical  and  experimental  investigation  of  the  oxidation  of  molybdenum  to 
molybdenum  trioxide  at  temperatures  at  which  the  trioxide  is  volatile  has  been  carried 
out.  A  rate  equation  has  been  derived  which  includes  the  effects  of  mass-transfer  param¬ 
eters  as  well  as  chemical  parameters  and  is  applicable  to  any  situation  for  which  a  mass- 
transfer  coefficient  can  be  obtained.  The  derived  equation  has  been  applied  to  previously 
reported  rate  data  taken  under  laboratory  conditions  using  air  and  oxygen-helium  mix¬ 
tures  flowing  at  low  subsonic  velocities,  and  a  set  of  values  for  the  activation  energy  of 
the  oxidation  of  molybdenum  has  been  obtained.  The  average  of  these  values  is 
8.51  x  107  J-(kg-mole)"1  (20.3  kcal-(g-mole)-1)  with  a  standard  deviation  of  2.1  percent. 

I  In  addition/32  molybdenum  hemisphere-cylinder  specimens  have  been  oxidized  in  a 
I  If  O  ) "  f  9 

Mach  2.1  heated  airstream  at  a  free-stream  stagnation  pressure  of  1.07  X  10°  N-m“ 

(10.5  atm)  and  free-stream  stagnation  temperatures  from  1765°  K  to  2120°  K.  The 
stagnation-point  surface  temperature  and  recession  rate  were  determined  for  each  of  the 
specimens.  The  observed  recession  rates  were  compared  with  the  corresponding  values 
predicted  theoretically  by  the  derived  rate  equation  and  agreement  within  a  factor  of 
about  2  was  obtained  for  most  specimens.  However,  a  few  specimens  ignited,  exhibiting 
recession  rates  much  greater  than  those  predicted  theoretically.  The  high  recession 
rates  associated  with  ignition  were  found  to  be  the  result  of  a  sharp  rise  in  surface  tem¬ 
perature  which  caused  some  melting  of  the  specimens  in  addition  to  their  normal 
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INTRODUCTION 


Molybdenum  differs  from  most  other  metals  in  that  over  an  extended  temperature 
range  below  the  melting  point  of  the  metal  its  most  stable  oxide,  molybdenum  trioxide,  is 
gaseous  and,  hence,  nonprotective  (ref.  1).  The  melting  point  of  molybdenum  is  2890°  K, 
but  the  normal  boiling  point  of  its  trioxide  is  1428°  K.  Furthermore,  under  appropriate 
conditions  the  trioxide  may  be  almost  completely  volatilized  as  much  as  several  hundred 


degrees  below  its  normal  boiling  point.  Whatever  the  conditions,  there  is  a  large  tem¬ 
perature  range  over  which  a  molybdenum  surface  situated  in  an  oxygen  containing  envi¬ 
ronment  is  bare  and,  therefore,  vulnerable  to  direct  attack  by  oxygen. 

The  unprotected  condition  of  molybdenum  at  high  temperatures  leads  to  very  high 
oxidation  rates  (refs.  2,  3,  and  4).  These  high  oxidation  rates  impose  a  drain  on  the 
oxygen  supply  just  adjacent  to  the  molybdenum  surface  and,  consequently,  parameters 
affecting  the  replenishment  of  this  supply  -  mass-transfer  parameters  -  may  be 
expected  to  influence  the  observed  oxidation  rate.  It  has  been  shown  both  theoretically 
(refs.  2  and  5  and  Discussion  of  ref.  3)  and  experimentally  (refs.  2,  3,  4,  and  6)  that  this 
is  in  fact  the  case:  the  oxidation  rate  of  molybdenum  is,  in  general,  dependent  upon 
mass-transfer  parameters  as  well  as  chemical  parameters,  although  under  certain  envi¬ 
ronmental  conditions  one  or  the  other  set  of  parameters  may  predominate. 

Because  the  oxidation  rate  of  molybdenum  at  high  temperatures  is  dependent  upon 
both  transport  and  chemical  parameters,  any  equation  which  attempts  to  describe  it  must 
include  the  effects  of  both  sets  of  parameters.  Such  an  equation  was  derived  theoreti¬ 
cally  in  reference  2  for  the  special  case  of  a  specimen  oxidizing  in  a  flowing  gas  stream 
under  conditions  such  that  laminar  boundary  layers  are  established.  The  effects  of  mass 
transfer  were  introduced  by  means  of  certain  simplifying  assumptions  which  resulted  in 
the  appearance  in  the  derived  rate  equation  of  an  empirical  proportionality  "constant" 
which  is  actually  a  function  of  various  environmental  parameters.  Thus,  although  the 
equation  derived  in  reference  2  serves  to  illustrate  the  interaction  of  transport  and  chem¬ 
ical  parameters  and  is  suitable  for  the  correlation  of  data  taken  under  similar  conditions, 
it  is  not  suitable  for  the  prediction  of  oxidation  rates. 

A  theoretical  rate  equation  similar  to,  but  more  general  than,  that  presented  in  ref¬ 
erence  2  is  derived  in  the  present  report.  In  this  treatment  the  experimental  system  is 
not  specified  and  the  effects  of  mass  transfer  are  introduced  by  means  of  a  mass-transfer 
coefficient  rather  than  by  recourse  to  the  simplifying  assumptions  utilized  in  reference  2. 
The  new  rate  equation  retains  the  capabilities  of  that  presented  in  reference  2  and,  in 
addition,  makes  possible  the  prediction  of  oxidation  rates  for  conditions  for  which  the 
value  of  the  mass-transfer  coefficient  can  be  obtained.  The  activation  energy  for  the 
oxidation  of  molybdenum  must,  of  course,  also  be  known.  The  determination  of  both  the 
activation  energy  and  the  mass-transfer  coefficient  is  discussed  in  this  report. 

The  present  investigation  has  also  included  an  experimental  study  for  the  purpose  of 
obtaining  molybdenum  oxidation  rate  data  under  conditions  more  severe  than  the  static 
and  low  subsonic  conditions  reported  in  the  literature  to  date  (refs.  2,  3,  4,  and  6).  This 
study  has  involved  the  oxidation  of  molybdenum  hemisphere- cylinders  in  heated  air- 
streams  at  a  Mach  number  of  2.1  and  a  free- stream  stagnation  pressure  of 
1.07  X  106  N~m~2  (10.5  atm).  The  average  steady-state  stagnation-point  surface 


2 


temperature  of  the  test  specimens  varied  from  1535°  K  to  2405°  K.  The  experimentally 
determined  oxidation  rates  were  compared  with  those  predicted  theoretically  by  the  rate 
equation  derived  in  this  report. 

SYMBOLS 

Wherever  necessary  the  symbols  and  units  used  in  the  various  references  cited  in 
this  report  have  been  changed  to  conform  with  those  given  in  this  section. 

a  speed  of  sound,  meters-second"! 

ao,al’a2  Parameters  in  equation  (B6) 

_  3 

C  concentration  of  oxygen,  kilograms- meter 

Cp  frozen  specific  heat  at  constant  pressure,  joules-kilogram"l-0K_1 

D  diffusivity,  meters^- second"  * 

D0  diffusivity  coefficient  (see  eq.  (B3)) 

d  diameter,  meters 

Ea  activation  energy,  joules- (kilogram- mole)- * 

H  enthalpy,  joules-kilogram"1 

h  heat-transfer  coefficient,  kilograms-meter"  -second" 

hj)  mass-transfer  coefficient,  kilograms- meter"  second"  1 

K  geometric  index  (see  eq.  (Al)) 

k  oxidation  rate  at  a  point  on  a  surface  (expressed  as  oxygen  consumed), 

kilograms-meter"  2- second"-^ 

kav  oxidation  rate  averaged  over  a  surface  (expressed  as  oxygen  consumed), 

kilograms-meter"  ^-second"! 
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length  of  specimen,  meters 
recession  rate,  meters- second"  1 
molecular  weight 

-  2  1 

mass-transfer  rate,  kilograms-meter  -second" 


NLe 

Lewis  number,  - 

NSc 

%a 

Mach  number 

nNu 

hcDx 

Nusselt  number,  — £— 

K 

Npr 

/iCp 

Prandtl  number,  — — 

K 

NRe 

puex 

Reynolds  number,  — — 

Nr6  w  modified  Reynolds  number, 


Pwuex 


Schmidt  number, 


temperature  exponent  (see  eq.  (B3)) 

mass-transfer  parameter,  assumed  in  reference  2  to  be  constant 

p 

pressure,  newtons-meter 

standard  pressure,  newtons- meter" 2  (see  eq.  (B3)) 
aerodynamic  heat  flux,  joules- meter" ^-second" 
universal  gas  constant,  joules- (kilogram-mole)"  1-°K"1 
radius  of  cross  section  of  body  of  revolution,  meters 
temperature,  °K 


standard  temperature,  °K  (see  eq.  (B3)) 


T  average  surface  temperature,  °K 

^30-40  average  of  stagnation-point  surface  temperature  at  30  and  40  seconds,  °K 
U  average  speed  of  oxygen  molecules,  meters- second"  * 

u  component  of  velocity  parallel  to  surface,  meter s- second"  * 

V  molar  volume  of  a  gas  at  its  normal  boiling  point,  meters^- (kilogram- mole) 

v  component  of  velocity  normal  to  surface,  meters- second-! 

w  mass  fraction 

X  mole  fraction 

x  coordinate  parallel  to  surface,  meters 

y  coordinate  normal  to  surface,  meters 

6q  concentration  boundary- layer  thickness,  meters 

6  normalized  enthalpy,  defined  by  equation  (A5) 

k  thermal  conductivity,  joules- meter-!- second"  !-°K-^ 

p,  viscosity,  kilograms-meter-!-second-! 

O 

p  density,  kilograms-meter-0 

cp  normalized  mass  fraction,  defined  by  equation  (A6) 

Subscripts: 

e  condition  at  edge  of  boundary  layer  or  other  mass-transfer  path 

He  helium 


x 


numerical  index 


Mo  molybdenum 

O2  oxygen 

t  stagnation  condition 

w  condition  adjacent  to  wall 

THEORETICAL  ANALYSIS 
Derivation  of  Rate  Equation 


It  has  been  shown  by  Gulbransen,  Andrew,  and  Brassart  (ref.  3)  that,  at  tempera¬ 
tures  at  which  the  reaction  products  are  volatile,  the  rate  determining  mechanism  for  the 
oxidation  of  molybdenum  involves  mobile  adsorption  of  oxygen  on  the  molybdenum  sur¬ 
face.  The  rate  equation  for  this  mechanism  according  to  Laidler,  Glasstone,  and  Eyring 
(ref.  7)  is 


By  making  use  of  the  expression  given  by  the  kinetic  theory  of  gases  for  the  average 
molecular  speed  U  (ref.  8),  this  equation  can  be  rewritten  as 


(2) 


In  this  form  the  rate  equation  for  mobile  adsorption  can  be  seen  to  be  identical  to  that 
derivable  from  simple  collision  theory.  (See  Discussion  of  ref.  3.) 

An  alternate  expression  for  equation  (2)  in  terms  of  the  mass  fraction,  rather  than 
the  concentration,  of  oxygen  at  the  surface  is 


w02,wPwu 
k  = - - -  ex] 


(3) 


For  steady- state  conditions  the  rate  of  consumption  of  oxygen  at  the  surface  must 
equal  the  rate  of  transport  of  oxygen  to  the  surface  -  that  is, 


(4) 
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The  oxygen  transport  rate  may  be  expressed  in  terms  of  a  mass-transfer  coefficient  as 


™02,w  =  hD(w02,e  -  w02,w) 


Therefore, 

k  =  hD(w02,e  "  w02,w) 


(5) 

(6) 


The  parameter 

unknown.  Therefore, 
tion  is 


w0  which  appears  in  equations  (3)  and  (6)  will,  in  general,  be 
it  is  eliminated  by  combining  the  two  equations;  the  resulting  equa- 


(7) 


Equation  (7)  is  a  general  expression  for  the  rate  of  oxidation  of  molybdenum  at  tem¬ 
peratures  at  which  its  trioxide  is  volatilized;  it  includes  the  effects  of  both  chemical  and 
transport  parameters.  Under  the  appropriate  conditions  either  of  these  two  sets  of 
parameters  may  predominate;  for  such  conditions  one  or  the  other  of  two  limiting  modi¬ 
fications  of  equation  (7)  will  result.  If  the  transport  term  hj y/pw  is  much  larger  than 
U  (  -Ea\ 

the  chemical  term  expl  j>  equation  (7)  reduces  to 

k  -  w02,efw  f 

which  is  the  limiting  equation  for  chemical  control. 

On  the  other  hand,  if  the  transport  term  is  much  smaller  than  the  chemical  term, 
equation  (7)  becomes 


k  =  w02,ehD  (®) 

which  is  the  limiting  equation  for  transport  control. 

In  the  application  of  equations  (7)  and  (9)  to  specific  problems  it  is  extremely  dif¬ 
ficult  to  obtain  creditable  values  for  the  parameters  pw  and  h^  if  the  pressure  and 
mean  molecular  weight  of  the  gas  adjacent  to  the  molybdenum  surface  are  significantly 
different  from  their  values  in  the  free  stream  or  reservoir.  (The  designation  of  the 
undepleted  main  body  of  the  gas  as  the  "free  stream"  or  the  "reservoir"  depends  upon 
whether  the  gas  is  flowing  or  essentially  static.)  The  conversion  of  oxygen  to  molybdenum 
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trioxide  by  reaction  with  the  molybdenum  surface  inevitably  results  in  changes  in  the 
pressure  and  mean  molecular  weight  of  the  gas  mixture  adjacent  to  the  surface.  These 
changes  are  negligibly  small  if  the  reaction  rate  is  chemically  controlled.  For  more 
rapid  reaction  rates  the  extent  of  the  changes  can  be  reduced  by  the  incorporation  of 
inert  diluents  in  the  oxidizing  gas  mixture.  If  a  sufficiently  large  fraction  of  inert  con¬ 
stituents  is  used,  the  pressure  and  mean  molecular  weight  of  the  gas  mixture  can  be 
maintained  effectively  uniform  even  when  the  reaction  proceeds  rapidly  enough  to  be  pre¬ 
dominantly  transport  controlled.  Because  of  these  considerations  it  will  be  assumed 
throughout  this  report  that  the  pressure  and  mean  molecular  weight  of  the  oxidizing  gas 
mixture  are  uniform  in  all  cases. 

In  view  of  the  foregoing  assumption,  equation  (7)  can  be  rewritten  for  the  special 
case  of  systems  where  the  temperature  of  the  gas  mixture  is  uniform.  For  such  systems 

Pe  =  Pw  =  P  (10) 

and  equation  (7)  becomes 


k  = 


-  —  exp 

'e  p  4  expi 


(ID 


For  systems  for  which  the  temperature  is  not  uniform  the  more  general  equation  (7)  must 
be  used. 

Equation  (11)  is  very  similar  to  a  rate  equation  derived  in  reference  2  for  the  spe¬ 
cial  case  of  oxidation  in  a  flowing  gas  stream  of  uniform  density  under  conditions  such 
that  laminar  boundary  layers  are  established.  The  two  equations  differ  only  in  the  trans¬ 
port  term  which  is  hj -^/p  in  equation  (11)  and  is  in  the  equation  derived  in 

reference  2.  It  can  be  seen  that  the  parameter  Pj  which  is  assumed  in  reference  2  to 
be  a  constant  is  actually  a  variable  defined  by  the  equation 


P 


1  “ 


hP6C 

pD 


(12) 


All  the  parameters  in  equations  (7)  and  (11),  with  the  exception  of  the  activation 
energy,  can  be  calculated  for  many  situations  of  interest.  (See  appendixes  A  and  B.) 
Therefore,  once  the  activation  energy  has  been  determined,  these  equations  can  be  used 
for  the  prediction  of  oxidation  rates  for  such  situations.  The  determination  of  the  activa¬ 
tion  energy  from  previously  reported  data  is  discussed  in  the  next  section. 
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Determination  of  Activation  Energy  From  Average  Oxidation- Rate  Data 

for  Flat- Plate  Specimens 

Data  on  the  oxidation  rate  of  molybdenum  taken  under  controlled  laboratory  condi¬ 
tions  are  reported  in  reference  2.  The  specimens  tested  were  thin  flat  plates  of  molyb¬ 
denum  which  were  oxidized  in  a  tube  furnace  in  slowly  flowing  gas  streams  with  flow 
velocities  of  1.76  to  8.00  m-sec'1.  The  specimens  were  all  alined  parallel  to  the  flow. 
The  exhaust  pressure  of  the  system  was  atmospheric  for  all  tests.  The  surface  tempera¬ 
ture  of  the  specimens  ranged  from  1069°  K  to  1644°  K.  Three  different  gas  mixtures 
were  used:  Air;  21.5  percent  O2,  78.5  percent  He;  and  13.6  percent  O2,  86.4  percent  He. 

The  reaction  conditions  and  gas  mixtures  employed  in  the  study  reported  in  refer¬ 
ence  2  were  such  that  the  density  could  be  considered  uniform  throughout  the  boundary 
layer  and,  thus,  the  appropriate  rate  equation  is  equation  (11).  However,  since  the  oxida¬ 
tion  rates  reported  in  reference  2  are  averages  over  the  entire  surface  of  the  specimens, 
rather  than  values  for  particular  points  on  the  surface,  equation  (11)  must  be  modified 
somewhat  in  order  to  be  applicable  to  these  data.  This  modification  is  accomplished  in 
the  following  paragraphs. 

Since  the  oxidation  rate  is  very  nearly  constant  over  any  narrow  strip  across  the 
surface  of  a  specimen  which  is  located  a  distance  x  from  the  leading  edge,  the  mean 
value  theorem  of  the  integral  calculus  gives  the  average  oxidation  rate  as 

kav  -  f  §*  k  (13) 

To  utilize  equation  (13)  it  is  necessary  to  have  an  expression  for  k  as  a  function 
of  x  which,  in  turn,  requires  an  expression  relating  hD  and  x.  In  appendix  A,  hD 
for  the  present  case  is  shown  to  be  given  by  the  equation 


(14) 


From  the  definition  of  the  Reynolds  number  and  equations  (11),  (13),  and  (14), 
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Integration  of  equation  (15)  produces  the  desired  expression  for  the  average  oxida¬ 
tion  rate: 


kav  - 


2(0.332)DwCe  J  0.332DW 


J£L\ 

pue) 


1/2 


*- 


M_\ 

puej 


1/2  77 


log. 


U 


-E, 


exp 


1  + 


fjjj_\1/2  U 
\^ue 


-E, 


exp 


RT 


Wy 


0.332DW  s3/N^q 


}  (16) 


RT 


w/ 


J 


Application  of  equation  (16)  to  the  data  presented  in  reference  2  yields  a  set  of 
values  for  the  activation  energy,  one  value  for  each  test.  (See  table  II.)  Only  the  data 
for  temperatures  of  1300°  K  and  above  have  been  considered,  since  these  are  the  data  for 
which  the  molybdenum  trioxide  reaction  product  was  completely  volatilized. 


EXPERIMENTAL  STUDY 


Facility  and  Test  Conditions 

All  tests  performed  in  this  study  were  carried  out  in  the  Langley  11-inch  ceramic- 
heated  tunnel  (ref.  9).  For  all  tests  the  oxidizing  gas  was  air,  the  stagnation  pressure  of 
the  airstream  was  1.07  x  106  N-m"2  (10.5  atm),  and  the  Mach  number  was  2.1.  The  free- 
stream  stagnation  temperature  varied  among  the  tests  from  1765°  K  to  2120°  K;  the 
values  for  the  individual  tests  are  presented  in  table  I. 

Test  Specimens 

All  specimens  tested  were  solid  0.0095-m-diameter  hemisphere- cylinders  which 
were  machined  from  sintered  molybdenum  of  a  purity  of  99.98  percent.  The  specimens 
were  mounted  in  a  water-cooled  sting  which  could  be  inserted  into  and  retracted  from 
the  heated  airstream  while  the  tunnel  was  in  operation. 

Photographs  of  a  typical  specimen  before  testing  are  presented  in  figure  1.  Fig¬ 
ure  1(a)  shows  the  specimen  unmounted  and  figure  1(b)  shows  the  specimen  mounted  in 
the  detachable  copper  nose  of  the  water-cooled  sting. 

Temperature  Measurement 

The  surface  temperature  of  the  specimens  tested  in  this  investigation  was  measured 
with  a  photographic  pyrometer.  The  theory  and  technique  of  photographic  pyrometry  are 
discussed  in  references  10  and  11.  The  pyrometer  used  in  the  present  study  was  the 
prototype  model  described  in  reference  10.  In  this  reference  the  precision  of  this  instru¬ 
ment  is  estimated  to  be  about  ±2  percent.  In  practice  it  has  been  found  that  precision  of 
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this  order  is  attainable  only  under  idealized 
conditions.  A  more  realistic  estimate  of  the 
precision  of  the  surface  temperatures  pre¬ 
sented  in  this  report  (table  I)  is  ±4  to 
±5  percent. 

Experimental  and  Analytical  Techniques 

The  operational  procedure  for  most  of 
the  tests  made  in  this  investigation  consisted 
of  starting  the  tunnel,  allowing  time  for  stable 
operation  to  be  achieved,  and  then  rapidly 
inserting  a  test  specimen  into  the  heated  air- 
stream  where  it  was  rigidly  held  for  a  pre¬ 
determined  length  of  time  at  the  end  of  which  it 
was  rapidly  retracted.  The  only  exceptions  to 
this  procedure  involved  the  retraction  before 
the  completion  of  their  scheduled  test  time  of 
a  few  specimens  which  unexpectedly  ignited. 

For  each  specimen  tested  the  surface 
temperature  and  recession  rate  were  deter¬ 
mined.  Because  both  of  these  parameters 
varied  across  the  surface  of  the  specimens, 
they  were  determined  at  a  single  common 
point  -  the  stagnation  point. 

The  stagnation-point  surface  tempera¬ 
ture  of  each  specimen  was  measured  at  30  sec 
after  insertion  of  the  specimen  into  the  heated 
airstream  and  every  10  sec  thereafter  until  the 
end  of  the  test.  Preliminary  tests  in  which 
temperatures  were  read  at  frequent  intervals 
up  to  30  sec  showed  that  the  stagnation-point 
surface  temperature  rose  rapidly  from  the 
ambient  value  for  the  first  few  seconds  and 
then,  except  for  the  few  specimens  that  ignited, 
gradually  tapered  off  to  a  steady-state  value 
which  was  generally  reached  between  20  and 
30  sec.  The  temperature  then  remained  rela¬ 
tively  constant  (within  ±5  percent  for  most 
specimens)  until  the  end  of  the  test. 


0  cm  I 

lb)  Mounted.  L-64-11,723 

Figure  1.-  Typical  hemisphere-cylinder  test 
specimen. 
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The  recession  rates  of  the  specimens  were  not  measured  directly  but  were  deter¬ 
mined  from  their  total  recessions  by  a  technique  which  is  described  in  the  following  par¬ 
agraphs.  The  total  recession  of  each  specimen  was  determined  by  measuring  its  total 
length  with  a  micrometer  before  and  after  testing. 

At  the  steady- state  temperatures  involved  in  this  investigation  the  rate  of  oxidation 
of  molybdenum,  at  a  given  temperature,  is  linear  with  respect  to  time.  However,  because 
there  was  a  significant  transient  period  before  the  specimens  reached  their  steady- state 
temperatures,  meaningful  recession  rates  could  not  be  calculated  simply  by  dividing  the 
total  recession  of  each  specimen  by  its  test  time;  some  compensation  for  the  recession 
which  occurred  during  the  transient  period  was  necessary.  This  compensation  was 
achieved  by  running  a  set  of  seven  40- sec  tests  at  various  free- stream  stagnation  tem¬ 
peratures.  (See  table  1(a).)  (Although  all  the  specimens  previously  tested  reached  their 
steady-state  temperatures  within  30  sec,  an  additional  10-sec  margin  was  allowed.)  The 
average  of  the  stagnation- point  surface  temperatures  at  30  and  40  sec,  designated 
T30  4Q>  was  determined  for  each  specimen.  In  figure  2  the  total  stagnation-point  reces¬ 
sion  of  each  of  the  seven  specimens  is  plotted  against  the  corresponding  value  of  T30..4O. 
Although  some  scatter  is  evident  among  the  data,  a  linear  trend  is  apparent  and,  conse¬ 
quently,  a  straight  line  was  fitted  to  the  data  by  the  least  squares  method.  The  ordinate 
of  each  point  on  this  line  was  then  assumed  to  represent  the  recession  in  the  first  40  sec 
of  any  test  specimen  with  a  value  of  T2Q  4Q  equal  to  the  abscissa  of  the  point. 


Figure  2.-  Variation  of  transient  recession  with  average  stagnation-point  surface  temperature. 


Twenty-two  specimens  were  tested  for  periods  longer  than  40  sec;  twenty-one  of 
these  oxidized  normally  and  did  not  ignite.  For  each  of  these  specimens  T30_4o  and> 
also,  the  average  stagnation -point  surface  temperature  from  40  sec  to  the  end  of  the  test 
T  were  determined.  (For  brevity,  the  latter  temperature  hereinafter  will  be  referred 
to  simply  as  the  "average  surface  temperature.")  From  figure  2  the  transient  recession 
of  each  test  specimen  corresponding  to  its  T30„4o  value  was  obtained.  This  transient 
recession  was  then  subtracted  from  the  specimen's  total  recession  and  the  difference  was 
divided  by  the  total  test  time  minus  40  sec.  The  resultant  value  was  assumed  to  be  the 
steady-state  recession  rate  of  the  specimen  corresponding  to  the  average  surface  temper¬ 
ature  T. 

The  technique  just  described  is  admittedly  less  precise  than  are  the  various  tech¬ 
niques  conventionally  employed  in  the  laboratory  for  the  determination  of  reaction  rates, 
techniques  such  as  continuous  weighing,  resistance  change  measurement,  and  oxygen  loss 
monitoring.  However,  none  of  these  conventional  techniques  were  applicable  to  the 
environment  used. 

It  can  be  seen  from  figure  2  that  the  actual  recession  of  a  specimen  in  40  sec  may 
deviate  as  much  as  0.00017  m  from  the  value  given  by  the  least  squares  line.  From 
table  1(b)  and  figure  2  it  can  be  determined  that  the  difference  between  the  transient  reces¬ 
sion  given  by  the  least  squares  line  and  the  total  recession,  for  the  specimens  tested 
longer  than  40  sec,  varies  from  0.00035  m  to  0.00310  m.  Thus,  the  uncertainty  in  the 
recession  rates  determined  varies  among  the  specimens,  ranging  from  at  least  5  percent 
to  49  percent  or  greater. 

The  technique  used  to  determine  the  recession  rates  of  the  specimens  which  oxi¬ 
dized  normally  required  significant  modification  for  the  determination  of  recession  rates 
due  to  ignition.  The  rapid  deterioration  which  resulted  from  ignition  was  in  every 
instance  accompanied  by  a  rapid  rise  in  surface  temperature  to  a  value  more  than  500°  K 
above  the  average  preignition  temperature.  Thus,  the  entire  preignition  period,  which 
varied  in  duration  from  23  to  130  sec,  constituted  a  transient  period  for  which  compensa¬ 
tion  was  necessary. 

This  compensation  was  achieved  as  follows.  The  duration  and  average  temperature 
of  the  preignition  period  were  determined  from  the  pyrometric  film  records  for  each 
specimen  which  ignited.*  (The  temperature  was  averaged  by  numerical  integration  using 

*In  one  test  in  which  ignition  occurred  the  specimen  (specimen  29)  was  not  immedi¬ 
ately  retracted  but  was  allowed  to  remain  in  the  heated  airstream  for  observation.  The 
specimen  receded  rapidly  for  about  21  sec  after  the  onset  of  ignition  and  then  reverted 
back  to  normal  oxidation  until  the  end  of  the  test.  The  duration  and  average  temperature 
of  this  postignition  period  of  normal  oxidation  were  included  with  the  specimen's  cor¬ 
responding  preignition  values. 
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the  trapezoid  rule;  temperatures  below  1428°  K  were  not  considered.)  The  theoretical 
recession  rate  associated  with  each  specimen's  average  preignition  temperature  was 
obtained  from  equation  (7).  This  rate  was  then  multiplied  by  the  specimen's  preignition 
time  to  yield  the  theoretical  preignition  recession.  Each  specimen's  ignition  recession 
rate  was  then  determined  by  subtracting  its  theoretical  preignition  recession  from  its 
measured  total  recession  and  dividing  the  difference  by  the  duration  of  ignition.  The 
ignition  recession  rates  determined  by  this  technique  are  probably  low  for  the  following 
reasons:  (1)  The  nose  of  each  specimen  which  ignited  assumed  an  irregular  configura¬ 
tion.  In  no  case  was  the  center  of  the  nose  the  farthest  protruding  point  when  the  speci¬ 
men  lengths  were  measured  after  testing.  Thus,  at  least  one  point  of  each  specimen,  the 
center  of  its  nose,  underwent  a  greater  recession  than  the  specimen's  measured  total 
recession.  (2)  As  pointed  out  in  the  section  "Results  and  Discussion,"  equation  (7),  the 
rate  equation  from  which  each  specimen's  preignition  recession  was  theoretically  deter¬ 
mined,  yields  values  higher  than  those  which  have  been  observed  experimentally  in  the 
preignition  temperature  range.  (3)  The  pyrometric  film  records  from  which  the  ignition 
times  were  determined  yield  readings  about  1  sec  apart  and,  therefore,  the  duration  of 
ignition  is  uncertain  in  each  case  by  about  2  sec.  The  ignition  durations  reported  in 
table  1(c)  are  the  maximum  values  derivable  from  the  film  records. 

RESULTS  AND  DISCUSSION 
Determination  of  Activation  Energy 

Table  II  reproduces  from  reference  2  certain  data  regarding  the  oxidation  of  thin 
flat  plates  of  molybdenum  to  molybdenum  trioxide  at  temperatures  and  flow  rates  at  which 
the  trioxide  is  completely  volatilized.3  Also  presented  in  this  table  are  values  of  the 
activation  energy  for  this  reaction  calculated  from  these  data  by  means  of  equation  (16); 
one  value  of  the  activation  energy  is  presented  for  each  data  point.  The  data  presented 

^  Equation  (7)  as  presented  yields  the  oxidation  rate  in  kilograms  of  oxygen  con¬ 
sumed  per  meter^  of  molybdenum  surface  per  second.  This  oxidation  rate  was  converted 
to  the  recession  rate  of  the  surface  in  meters  per  second  by  the  following  equation: 

2mMo 

"  3mO/Mo 

O 

For  consistency  with  the  conventions  of  the  present  report  the  values  of  these 
data  have  been  converted  to  the  mksK  system  of  units  and  the  oxidation  rates  have  been 
expressed  in  terms  of  the  mass  rate  of  oxygen  consumed  rather  than  the  mass  rate  of 
molybdenum  lost. 


14 


in  table  II  are  divided  among  five  experimental  series  comprising  five  different  mass 
flow  rates  and  three  gas  mixtures:  air  and  two  oxygen-helium  mixtures.  As  demon¬ 
strated  in  reference  2  all  oxidation  rates  presented  were  to  some  extent  transport 
controlled. 

Table  III  presents  the  average  value  of  the  activation  energy  for  each  test  series 
as  calculated  by  equation  (16)  and  also  as  calculated  by  a  corresponding  integrated  rate 
equation  derived  in  reference  2.  In  addition,  table  III  presents  the  overall  average  value 
of  the  activation  energy  and  the  percent  standard  deviation  from  this  average  for  each  of 
the  two  rate  equations  considered.  The  percent  standard  deviation  obtained  by  using 
equation  (16),  2.1  percent,  is  an  improvement  over  the  value  of  5.7  percent  obtained  with 
the  integrated  rate  equation  presented  in  reference  2.  (The  value  of  3.6  percent  given  in 
reference  2  is  the  mean  deviation  and  not  the  standard  deviation.) 

It  can  be  seen  from  table  III  that  not  only  does  the  use  of  equation  (16)  yield  more 
consistent  values  for  the  activation  energy  when  all  the  experimental  data  are  considered 
together  but  it  also  significantly  reduces  the  differences  among  the  average  values  for 
the  individual  series.  The  maximum  difference  in  the  average  value  of  the  activation 
energy  between  two  series  obtained  by  using  equation  (16)  is  0.16  x  107  J-(kg-mole)  , 
whereas  with  the  equation  in  reference  2  it  is  0.86  x  107  J-(kg-mole) 

The  overall  mean  value  of  the  activation  energy  as  calculated  by  equation  (16)  is 
8.51  x  107  J-(kg-mole)_1  (or  20.3  kcal-(g-mole)"1)  and  as  calculated  by  the  equation  from 
reference  2  it  is  8.95  X  107  J-(kg-mole)-1  (or  21.3  kcal-(g-mole)-1).  The  former 
agrees  within  3  percent  and  the  latter  agrees  within  8  percent  with  the  value 
8.25  x  107  J-(kg-mole)"1  (l9.7  kcal-(g-mole)-1)  obtained  by  Gulbransen,  Andrew,  and 
Brassart  (ref.  3)  from  an  Arrhenius  plot  of  experimental  data  taken  by  them.  In  general, 
the  use  of  this  type  of  plot  for  the  determination  of  the  activation  energy  of  a  material 
which  forms  a  volatile  oxide  is  unreliable  due  to  the  possible  effects  of  transport  param¬ 
eters  on  the  measured  oxidation  rates  (refs.  2  and  5).  However,  when  the  data  plotted 
have  been  obtained  under  predominantly  chemically  controlled  conditions,  this  source  of 
error  is  eliminated.  Gulbransen,  Andrew,  and  Brassart  were  able  to  obtain  some  data 
under  chemically  controlled  conditions  and  only  these  data  were  considered  in  their 
determination  of  the  activation  energy. 

It  is  apparent  from  the  results  just  presented  that  equation  (16)  reduces  the  data  of 
reference  2  more  successfully  than  does  the  integrated  rate  equation  presented  in  that 
reference. 


Recession-Rate  Data  for  Hemisphere- Cylinder  Specimens 

The  recession  rates  and  corresponding  average  surface  temperatures  of  the  molyb¬ 
denum  hemisphere-cylinder  specimens  tested  in  the  present  investigation  —  both  those 
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Figure  3.-  Variation  of  recession  rate  with  surface  temperature, 


that  behaved  normally  and  those  that  ignited  -  are  listed  in  table  I  and  plotted  in  figure  3. 
Also  plotted  in  figure  3  is  a  theoretical  curve  of  recession  rate  as  a  function  of  surface 
temperature  calculated  by  means  of  equation  (7)  and  based  on  the  specimen  configuration 
and  test  conditions  utilized.  This  curve  extends  from  1428°  K,  the  normal  boiling  point 
of  molybdenum  trioxide,  to  2890°  K,  the  melting  point  of  pure  molybdenum. 

It  can  be  seen  that  the  nonignition  data  broadly  follow  the  theoretical  curve  over 
the  temperature  range  covered  by  the  data  although  the  data  points  all  lie  somewhat 
below  the  curve.  The  degree  of  scatter  exhibited  by  the  data  is  consistent  with  the  esti¬ 
mated  uncertainty  of  the  measured  temperatures  and  recession  rates.  Furthermore, 
none  of  the  measured  recession  rates  of  specimens  which  did  not  ignite,  except  that  at 
1835°  K,  deviate  from  the  corresponding  values  predicted  by  equation  (7)  by  a  factor 
greater  than  3.3.  On  the  average  the  agreement  is  somewhat  better  than  this  -  a  factor 
of  about  2.  This  agreement  is  quite  reasonable,  as  high-temperature  kinetics  studies  go, 
even  under  more  idealized  circumstances;  for  example,  Gulbransen,  Andrew,  and 
Brassart,  studying  the  oxidation  of  molybdenum  under  laboratory  conditions  and  con¬ 
sidering  only  the  special  case  of  chemical  control,  also  obtained  differences  between 
theory  and  experiment  by  a  factor  of  about  3  (ref.  12).  It  should  be  noted  that  in  the  pres¬ 
ent  investigation  purely  chemical  control  was  not  maintained;  the  degree  of  transport 
influence,  calculated  by  using  equations  (7),  (8),  and  (9),  was  11  percent  at  1535°  K  and 
29  percent  at  1965°  K. 

In  contrast  with  the  behavior  of  those  specimens  which  oxidized  normally,  the 
experimentally  determined  recession  rates  of  the  few  specimens  which  ignited  are  much 
higher  than  the  corresponding  theoretical  values  calculated  by  means  of  equation  (7). 
Furthermore,  the  experimentally  determined  values  are  almost  certainly  low,  as  pointed 
out  in  the  section  "Experimental  and  Analytical  Techniques." 

The  behavior  of  the  specimens  that  ignited  can  best  be  understood  by  considering 
the  nature  of  the  phenomenon  of  ignition.  According  to  Hill,  Adamson,  Foland,  and 
Bressette  (ref.  13)  ignition  occurs  when  the  rate  of  heat  input  to  an  oxidizing  surface  due 
to  its  own  exothermic  reaction  exceeds  the  rate  of  heat  removal  from  the  surface.  A 
characteristic  of  ignition  is  a  rapid  increase  in  the  temperature  of  the  surface  until  some 
mechanism  —  frequently  melting  -  comes  into  effect  to  bring  the  heat  input  and  removal 
into  balance  again.  The  rapid  rise  in  temperature  which  is  characteristic  of  ignition  was 
observed  in  this  investigation;  this  is  shown  in  figure  4  in  which  temperature  is  plotted 
against  time  for  three  of  the  four  specimens  that  ignited  and  also  for  one  specimen  that 
reached  a  high  steady- state  temperature  but  did  not  ignite.  (A  temperature-time  record 
was  not  obtained  for  one  specimen  that  ignited.)  The  zero  of  time  as  plotted  in  figure  4 
is  not  absolute  for  any  of  the  specimens  but  has  been  chosen  in  each  case  so  as  to  produce 
the  maximum  overlap  among  the  data  points  at  lower  temperatures.  Significant  overlap 
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Figure  4.-  Variation  of  stagnation-point  surface  temperature  with  time. 

is  evident  up  to  about  1950°  K  but  above 
this  temperature  the  two  sets  of  data 
begin  to  diverge.  As  time  increases, 
the  temperature  of  the  normally  oxi¬ 
dizing  specimen  gradually  levels  out; 
the  temperature  of  the  specimens  which 
ignited  continues  to  increase  steadily 
for  a  short  time  until  at  about  2100°  K  a 
sharp  rise  of  195°  K  to  265°  K  occurs 
in  approximately  1  second. 

That  melting  occurred  with  the 
specimens  which  ignited  in  this  inves¬ 
tigation  is  illustrated  in  figures  5 
and  6.  In  figure  5  molten  material  can 
clearly  be  seen  flowing  back  from  the 
nose  of  a  specimen  which  had  just 
ignited.  In  figure  6  the  specimens 
which  ignited  can  be  seen  to  have  under¬ 
gone  significant  alteration  of  their 
original  configuration,  whereas  the  Figure  5 .  Specjmen  undergoing  ignition.  L.65_9024 
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specimens  which  did  not  ignite  retained  substantially  their  original  shape  despite  having 
undergone  measurable  recession.  Also,  two  specimens  which  ignited  have  deposits  of 
molten  material  on  their  sides  which  apparently  flowed  back  from  the  region  where 
melting  occurred  and  solidified  after  the  specimens  were  removed  from  the  heated  air- 
stream.  The  other  two  specimens  were  virtually  destroyed  and  the  deposit  of  molten 
material  on  their  sides  is  difficult  to  observe. 

Inasmuch  as  the  removal  of  molybdenum  from  the  nose  of  those  specimens  which 
ignited  was  caused  by  melting  as  well  as  by  oxidation,  it  is  not  surprising  that  the  meas¬ 
ured  recession  rates  of  these  specimens  are  considerably  higher  than  the  rates  predicted 
by  a  mechanism  that  considers  oxidation  alone. 

\  I  it  is  not  known  why  the  temperature  at  which  melting  occurred  was  determined  to 
be  about  500°  K  below  the  normal  melting  point  of  molybdenum*  Errors  in  temperature 
measurement  may  account  for  some  of  the  difference  since  the  presence  of  a  molten  film 
undoubtedly  changes  somewhat  the  emissivity  of  the  surface  observed.]  However,  a 


temperature  error  of  about  20  percent  is  unlikely.  It  is  probable  that  the  formation  of  a 
solution  of  metal  and  oxide  occurs  at  the  specimen  surface,  which  lowers  the  melting 
point  below  that  of  pure  molybdenum. 

Also  uncertain  is  the  reason  that  some  specimens  which  ignited  did  so  under  condi¬ 
tions  apparently  no  more  severe  than  those  experienced  by  other  specimens  which  did  not 
ignite.  (See  table  I.)  It  is  known  from  unpublished  tests  carried  out  in  the  Langley  arc- 
heated  materials  jet  that  at  sufficiently  severe  conditions  (Mach  2,  2535°  K  to  2870°  K 
free -stream  stagnation  temperature,  and  1.07  x  106  N-rn"2  stagnation  pressure)  ignition 
occurs  consistently.  On  the  other  hand,  none  of  the  tests  reported  in  references  2,  3,  4, 
and  6,  all  of  which  involved  subsonic  or  static  conditions  and  gas  temperatures  of  less 
than  2100°  K,  resulted  in  ignition.  Apparently  the  conditions  employed  in  the  present 
investigation  lie  in  a  transition  region  between  consistent  normal  oxidation  and  consistent 
ignition.  In  this  region  small  differences  in  heat  input  and/or  heat  removal  may  deter¬ 
mine  whether  a  specimen  ignites  or  behaves  normally. 

CONCLUDING  REMARKS 

A  rate  equation  for  the  oxidation  of  molybdenum  to  its  trioxide,  applicable  at  tem¬ 
peratures  at  which  the  trioxide  is  volatile,  has  been  derived  theoretically;  this  equation 
takes  into  account  the  effects  of  both  chemical  and  mass-transfer  parameters.  The  equa¬ 
tion  is  more  general  than  a  similar  equation  presented  in  NASA  Technical  Note  D-222 
and  should,  in  fact,  be  applicable  to  any  situation  where  the  mass-transfer  coefficient 
of  oxygen  to  a  molybdenum  surface  can  be  determined.  The  equation  has  been  applied 
to  low  flow  velocity  data  contained  in  the  aforementioned  Technical  Note  and  values  for 
the  activation  energy  of  the  oxidation  of  molybdenum  to  its  trioxide  have  been  obtained. 
The  standard  deviation  among  the  values  is  2.1  percent  and  the  average  value  of 
8.51  x  10?  J- (kg- mole) "1  (20.3  kcal-(g-mole)_l)  is  within  3  percent  of  the  value 
8.25  x  10?  J-(kg-mole)"1  (l9.7  kcal-(g-mole)-1)  determined  by  Gulbransen,  Andrew,  and 
Brassart  (Journal  of  the  Electrochemical  Society,  September  1963)  by  a  conventional 
Arrhenius  plot  using  data  obtained  under  predominantly  chemically  controlled  conditions. 

The  theoretically  derived  rate  equation  has  also  been  applied  to  the  oxidation  of 

9 

molybdenum  hemisphere-cylinders  in  heated  Mach  2.1  airstreams  at  1.07  x  10u  N-m~* 
free- stream  stagnation  pressure.  For  specimens  with  stagnation- point  surface  tempera¬ 
tures  remaining  below  about  2100°  K,  the  recession  rates  predicted  theoretically  by  the 
derived  rate  equation  agreed,  on  the  average,  within  a  factor  of  about  2  with  the  rates 
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determined  experimentally.  Specimens  with  stagnation-point  surface  temperatures 
rising  above  2100°  K  ignited.  This  ignition  was  found  to  be  characterized  by  a  rapid  rise 
in  specimen  surface  temperature  until  melting  occurred  and,  as  might  be  expected, 
resulted  in  recession  rates  considerably  in  excess  of  those  accountable  by  oxidation 
alone. 

Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Langley  Station,  Hampton,  Va.,  November  10,  1965. 


21 


APPENDIX  A 


CALCULATION  OF  MASS- TRANSFER  COEFFICIENTS 


Analogy  Between  Mass  Transfer  and  Heat  Transfer 

The  phenomenon  of  heat  transfer  across  an  enthalpy  or  temperature  gradient  has 
been  investigated  extensively  and  heat-transfer  coefficients  are  available  for  a  large 
variety  of  experimental  situations.  Mass-transfer  coefficients,  on  the  other  hand,  are 
less  readily  available.  However,  for  certain  situations  mass-transfer  coefficients  can  be 
obtained  from  known  heat- transfer  coefficients  by  means  of  analogies  which  exist  between 
mass  transfer  and  heat  transfer  for  these  situations.  Two  forms  of  mass-heat  transfer 
analogy  which  were  utilized  in  this  investigation  are  derived  and  discussed  in  this 
appendix. 

Consider  a  binary,  compressible,  laminar  boundary  layer  in  which  no  homogeneous 
chemical  reactions  occur.  The  governing  equations,  as  developed  in  reference  14,  are 
as  follows: 

Continuity  equation 


3  [  K 


A| 


pvrK)  =  0 


(Al) 


where  K  =  0  for  a  two-dimensional  body  and  K  =  1  for  a  body  of  revolution. 
Momentum  equation 


,  9u  3u 


5st±U 

dx  8yf  Sy 


(A2) 


Energy  equation 


p«s+ts).4j*_bua 

\  9x  9y/  fy'fJpr  sy/  9y 


Ml  - 


3  u 


Npr/9y\2 


sy 


pDh- 


H, 


Swi 


Nl  e/V 


(A3) 


Continuity  equation  for  each  species 


P 


u 


d\V± 

Ox 


+  V 


OWi 


3  (  in 

9y\Nsc  9y/ 


=  0 


(A4) 


Consider  the  last  two  terms  on  the  right-hand  side  of  equation  (A3).  The  next  to 
last  term  represents  heat  generated  by  viscous  shear  stresses  in  the  boundary  layer,  and 
the  last  term  represents  heat  transfer  by  diffusion.  Note  that  when  these  two  terms  are 
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negligible,  equations  (A3)  and  (A4)  are  of  similar  form.  If  these  two  equations  (with  the 
dissipation  and  diffusion  terms  of  equation  (A3)  neglected)  are  expressed  in  terms  of  the 
dimensionless  variables 


,Ht,w  -  Ht,e/ 


(A5) 


0i 


( wi  -  wi,w  ' 
\wi,e  "  wi,w/ 


(A6) 


the  similarity  of  form  is  maintained  and,  in  addition,  the  boundary  conditions  on  the 
dependent  variables  are  the  same  -  that  is, 


and 


/  90  90\  9  \i  90\ 

PiU  *  +  V  ayj  =  9y\Npr  WJ 

(A7) 

1  90i  90i\  9/  u  90 i\ 

p\  ~9x  +  V  ~dyj  =  "¥\Ns^  ~9y j 

(A8) 

0  =  0i  =  0 

(y  =  o) 

(A9) 

II 

\ 

(y  -  °°) 

(A10) 

If  the  Lewis  number  is  unity  (Npr  =  NScj,  equations  (A7)  and  (A8)  are  not  just  sim¬ 
ilar  in  form  but  are  identical  and,  hence,  6  =  6^. 

Therefore, 


/  _  /90i\ 

\9y/w  \  /  w 


(All) 


From  Fourier's  law  of  conduction  and  the  definition  of  the  heat-transfer  coefficient 


\  *Vw  w 


(A12) 


and  from  Fick’s  first  law  of  diffusion  and  the  definition  of  the  mass-transfer  coefficient 


^i,w  ~  "Pw^w(wi,e  wi,w) 


90i\  _ 

*/w" 


^D(wi,e  wi,w 


(A13) 
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Hence, 


hD  _  hcp,w 

P\v  ^  \\r  Kw 


(A14) 


Equation  (A14)  is  a  statement  of  the  mass-heat  transfer  analogy  for  situations  in 
which  simultaneous  mass  and  heat  transfer  are  involved.  It  must  be  remembered  that 
this  analogy  is  rigorously  valid  only  when  N^e  =  1  (note  that  this  condition,  in  addition 
to  specifying  the  equality  of  the  Schmidt  and  Prandtl  numbers,  eliminates  the  diffusion 
term  in  equation  (A3))  and  there  is  no  appreciable  heat  generation  in  the  boundary  layer 
due  to  viscous  sheer  stresses.  The  latter  condition  is  met  whenever  the  local  flow 
velocity  is  small  compared  with  the  velocity  of  sound.  This  is  the  case  for  all  points  on 
a  body  situated  in  a  low  subsonic  flow  stream  but  for  transonic  and  supersonic  flow 
streams  it  is  the  case  only  near  the  stagnation  point.  The  requirement  that  NLe  =  1  is 
seldom  met  exactly  but  is  generally  assumed  to  be  satisfied  for  values  between  about  0.7 
and  1.5.  Air  over  a  large  range  of  temperatures  and  pressures  meets  this  less  stringent 
requirement  as  do  a  number  of  other  gas  mixtures  (ref.  14). 

A  somewhat  modified  form  of  mass-heat  transfer  analogy  can  be  derived  for  the 
special  case  of  mass  transfer  through  a  boundary  layer  within  which  no  heat  transfer 
takes  place  without  assuming  the  Lewis  number  to  be  unity.  The  governing  equations  for 
this  special  case  are  equations  (Al),  (A2),  and  (A8).  Because  there  is  no  heat  transfer  in 
the  boundary  layer,  equation  (A7)  is  eliminated. 


In  the  development  of  this  form  of  analogy  it  is  necessary  to  consider  another  spe¬ 
cial  case  -  namely,  the  transfer  of  heat  through  a  boundary  layer  within  which  no  mass 
transfer  or  viscous  dissipation  takes  place.  For  this  case  equations  (Al),  (A2),  and  (A7) 
are  the  governing  equations  and  equation  (A8)  is  eliminated.  The  two  special  cases  can 
be  seen  to  be  similar  since  they  have  equations  (Al)  and  (A2)  in  common  and  equa¬ 
tions  (A7)  and  (A8)  differ  only  in  their  dimensionless  tei’ms  4>,  Npr,  and  Ngc)-  If 
both  the  Schmidt  number  and  the  Prandtl  number  are  independent  of  x  and  y,  then  any 
solution  to  equation  (A7)  can  be  transformed  into  a  solution  to  equation  (A8)  by  replacing 
the  Prandtl  number,  wherever  it  appears,  by  the  Schmidt  number.  Obviously,  a  similar 

transformation  exists  between  (4^)  and  (— i)  and,  hence,  in  view  of  equations  (A12) 

\ty/w  V  ®y  /  w 


and  (A13),  between 


‘D 


Pw^' 


and 


he 


p,w 


This  transformation  form  of  mass-heat  transfer 


w 


'W 


analogy  is  equivalent  to  that  presented  by  Eckert  in  reference  15. 


As  stated  previously  this  form  of  analogy  cannot  be  applied  to  situations  in  which 
both  mass  transfer  and  heat  transfer  occur  simultaneously  to  a  significant  extent.  In 
addition,  since  the  Prandtl  and  Schmidt  numbers  have  been  assumed  to  be  independent  of 
x  and  y,  the  composition  of  the  gas  must  be  relatively  homogeneous  throughout  the 
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boundary  layer.  This  assumption  is,  of  course,  met  when  the  rate  of  reaction  of  oxygen 
at  the  surface  is  slow  compared  with  the  replenishment  rate  or,  in  other  words,  when  the 
reaction  is  essentially  chemically  controlled.  When  chemical  control  does  not  apply, 
homogeneity  can  still  be  approximated  by  the  presence  of  a  large  percentage  of  inert  gas. 
Finally,  as  was  the  case  with  the  first  form  of  analogy  derived,  the  local  Mach  number 
must  be  relatively  small. 


Thin  Flat  Plates 


The  thin  flat-plate  specimens,  the  oxidation  rates  of  which  are  presented  in  refer¬ 
ence  2  and  reproduced  herein  in  table  H,  were  oxidized  under  conditions  such  that  heat 
transfer  within  the  boundary  layer  was  negligible.  The  specimens  were  oxidized  in 
flowing  gas  streams  with  very  low  flow  velocities,  of  the  order  of  1  to  10  m-sec  .  The 
gas  mixtures  used  were  air  and  two  helium-oxygen  mixtures,  all  of  which  had  at  least 
79  mole  percent  inert  constituents.  In  view  of  these  conditions,  mass-transfer  coeffi¬ 
cients  for  these  specimens  can  be  obtained  from  related  heat-transfer  coefficients  by 
means  of  the  second  form  of  mass-heat  transfer  analogy  derived. 


Heat-transfer  coefficients  (as  defined  in  this  report)  are  related  to  Nusselt  numbers 
by  the  following  equation: 


N 


Nu 


_  h^p,wx 


aw 


(A15) 


The  Nusselt  number  for  a  thin  flat  plate  in  laminar  flow  is  given  in  reference  15  as 


NNu  =  0.332  \^/Npr  (A16) 

Combining  equations  (A15)  and  (A16)  and  dividing  by  x  yields 

^  =  0432,3/5^^  (A17) 

According  to  the  applicable  analogy  equation  (A17)  may  be  transformed  as  follows: 

J±D_  =  \3/nZ  \/%^ 

pwDw  x  V  Sc  V  Re 

Therefore, 

hD  =  q.332  ^1^/N^v^ 


(A18) 


(A19) 
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Hemisphere-  Cylinders 

The  hemisphere- cylinders  tested  in  this  investigation  were  oxidized  under  condi¬ 
tions  such  that  the  amount  of  heat  transfer  which  accompanied  the  mass  transfer  was  not 
negligible.  The  only  gas  mixture  used  was  air  for  which  the  Lewis  number  at  the  temper¬ 
atures  and  pressures  involved  is  reasonably  close  to  unity  (ref.  14).  Oxidation  rates 
were  determined  only  at  the  stagnation  point  of  the  specimens.  These  considerations 
point  to  the  use  of  the  first  form  of  mass-heat  transfer  analogy  derived,  as  expressed  by 
equation  (A14). 

From  the  definition  of  the  Nusselt  and  Reynolds  numbers 

h^p,w 

Kw 

Combining  equations  (A14)  and  (A20)  gives 


^Nu.w  rMa,wVw 


1/2 


w 


Mwx 


(A20) 


hD  _  %u,w  |NMa,wawPw 
PwDw  ■  Vn^V  w 

In  reference  16  is  presented  a  plot  of  N^a  w  against  x/d  for  each  of  several 
free- stream  Mach  numbers  greater  than  1.  In  the  stagnation  region  all  the  curves 
asymptotically  approach  a  straight  line  passing  through  the  origin  with  a  slope  of  2.4. 
Therefore,  in  the  stagnation  region  for  free- stream  Mach  numbers  greater  than  1  the 
local  Mach  number  is  independent  of  the  free- stream  Mach  number  and  is  given  by  the 
expression 


a/2 


(A21) 


2.4x 


%a,w  _  (j 

Substitution  of  equation  (A22)  into  equation  (A21)  yields 


Therefore, 


hD  _  NNu,w  (2Aavjpw^2 
£\v^W  \/Nr6  w  \  ^W*!  / 


hD  = 


nNu,wDw| 

2-4aw^w3 

V 

^Re,w 

V  Mw^ 

(A22) 


(A23) 


(A24) 
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Reference  17  presents  a  series  of  plots  of  %[U)W/(/NRe)W  at  the  stagnation  point 
of  an  axially  symmetric  body  as  a  function  of  Tw/Te  for  each  of  several  values  of  the 
Prandtl  number.  These  plots  are  represented  quite  well  by  the  equation 


Therefore, 


NNu,w 

Re,w 


\0.9 


0.7  +  0.064 
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(A25) 
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(A26) 


Since  this  expression  for  the  mass-transfer  coefficient  is  relatively  insensitive  to 
t  /x  this  ratio  was  assumed  to  be  unity  in  the  determination  of  the  theoretical  curve 

Vvf 

presented  in  figure  3  of  this  report. 
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VALUES  OF  PARAMETERS  EMPLOYED  IN  COMPUTATIONS 

In  the  application  of  various  equations  presented  in  this  report  it  is  necessary  to 
have  values  for  the  parameters  Ea,  p,  U,  D,  p,  a,  Cp,  and  k.  The  values  and 
equations  for  these  parameters  which  were  employed  in  this  investigation  are  given  in 
this  appendix. 


Activation  Energy 

The  activation  energy  was  assumed  to  be  8.55  X  10^  J-(kg-mole)-^ 

(20.3  kcal-(g-mole)-l).  This  is  the  average  value  determined  by  means  of  equation  (16) 
from  the  data  presented  in  reference  2. 


Density 


The  equation  for  the  density  was  obtained  by  a  rearrangement  of  the  ideal  gas  equa¬ 
tion  —  that  is 


(Bl) 


The  molecular  weights  of  the  mixtures  were  taken  to  be  the  averages  of  the  molecular 
weights  of  the  individual  gases,  weighted  according  to  their  mole  fractions. 


Average  Velocity  of  Oxygen  Molecules 

The  equation  used  for  the  average  velocity  of  oxygen  molecules  is  that  derived  in 
the  kinetic  theory  of  gases.  (See,  for  example,  ref.  8.)  It  is 


U  = 


'  8RT  \1//2 

?mo2) 


(B2) 


Diffusivity 

A  semiempirical  equation  for  the  diffusivity  of  various  binary  gas  mixtures  is 
given  in  reference  18  -  namely, 


D  =  D0 


Po 

P 


(B3) 


The  standard  temperature  T0  and  standard  pressure  pQ  employed  in  reference  18 
are  273°  K  and  1.013  X  105  N  -m~2  (1  atm),  respectively.  The  constants  D0  and  n  are 
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given  for  a  number  of  gas  mixtures,  including  air.  Substitution  of  the  values  given  for 
T0,  p0,  D0,  and  n  into  equation  (B3)  and  multiplication  by  a  conversion  factor  to 

yield  units  of  meters^  per  second  produces  the  equation 


D 


1.00  X  10-4T1-75 


air 


(B4) 


Reference  18  does  not  give  the  values  of  D0  and  n  for  oxygen-helium  mixtures.  How¬ 
ever ,  n  was  assumed  to  be  1.75  since  this  is  the  value  given  for  nearly  two-thirds  of 
the  mixtures  cited,  including  the  only  mixture  containing  helium  (helium-argon).  A  value 
for  D0  was  obtained  for  the  oxygen-helium  mixtures  by  means  of  a  relationship  pre¬ 
sented  in  reference  19.  In  this  reference  the  diffusivity  of  binary  gas  mixtures  is  repre¬ 
sented  by  an  expression  similar  to  equation  (B3).  The  term  corresponding  to  D0  is 


proportional  to  the  parameter 


+  V, 


1/3 


Evaluation  of  the  ratio 


of  the  value  of  this  parameter  for  one  gas  mixture  to  its  value  for  a  second  gas  mixture 
should  yield  the  ratio  of  the  values  of  D0  for  the  two  mixtures  and,  hence,  if  the  value 
of  D0  is  known  for  one  of  the  mixtures,  its  value  for  the  other  can  be  determined. 
Employing  this  scheme  and  using  the  known  value  of  D0  for  air  gives  the  following 
equation: 


D 


02,He 


2.22  x  10 


-4  T 


1.75 


(B5) 


Viscosity 

Reference  18  gives  the  following  semiempirical  equation  for  the  viscosity  of  certain 
gases: 


M  = 


a0T 


1/2 


1  + 


(B6) 


If  the  values  for  the  constants  a0,  aj,  and  a2  which  are  given  in  this  reference  for  air 
are  substituted  into  equation  (B6)  and  an  appropriate  factor  is  introduced  to  convert  the 
viscosity  from  poise  to  newtons  per  meter^  per  second,  the  following  equation  is  obtained: 


^air  “ 


1.488  X  10~6T1/2 


1  + 


122.1  X  10 


(B7) 


Values  of  the  constants  in  equation  (B6)  are  not  given  in  reference  18  for  oxygen- 
helium  mixtures  and,  hence,  this  equation  could  not  be  applied  to  these  mixtures. 
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However,  reference  19  gives  a  general  equation  for  the  viscosity  of  a  binary  gas  mixture 
which  is  applicable,  in  principle,  to  any  pair  of  gases  with  known  individual  viscosities. 
This  equation  is 


ll12~ 


H 


1  + 


1  X2 


l  +  (qi/M2)1/2(m2/ml)l/4 


2\/2  Xi  ;  *1/2 

(1  +  mj/mg) 


V2 


1  + 


1  *1 


[l  +  (M2/Mi)l/2(m1/m2)l/^ 


2V2  x2 


(1  +  m2/m1) 


1/2 


(B8) 


The  individual  viscosities  of  oxygen  and  helium  were  not  obtained  from  equa¬ 
tion  (B6)  although  the  values  of  a0,  a1?  and  a2  for  these  gases  are  presented  in  ref¬ 

erence  18,  because  the  temperature  range  of  applicability  cited  for  these  two  gases  is 
too  restricted. 

The  equation  used  for  the  viscosity  of  helium  was  obtained  from  reference  20. 
With  appropriate  conversion  of  units  it  is 


^He  “ 


4.230  x  10~7T3//2 
t0.826  _  0.409 


(B9) 


Sutherland's  equation  (see  ref.  19)  was  used  for  the  viscosity  of  oxygen.  The  con¬ 
stants  in  Sutherland's  equation  were  obtained  by  fitting  the  equation  to  the  viscosity  values 
presented  in  reference  21.  The  resulting  equation  is 


_  1.612  x  iq-St1/2 
M°2  "  17.95 

X  *f*  rp 


(BIO) 


Speed  of  Sound 

The  equation  used  for  the  speed  of  sound  in  air  is 

aair  =  25.30TW0‘4636  (Bll) 

This  equation  is  an  empirical  fit  to  tabulated  values  presented  in  reference  21. 
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Specific  Heat 

The  equation  used  for  the  specific  heat  of  air  at  constant  pressure  is 

5p,air  =  9.781  x  1C)2  +  0.1718TW  (B12) 

Equation  (B12)  is  also  an  empirical  fit  to  tabulated  values  presented  in  reference  21. 


Thermal  Conductivity 

The  equation  used  for  the  thermal  conductivity  of  air  is 

*air  =  8>476  X  10"3  +  6>266  X  10_5tw  '  8.501  x  10"9Tw2  (B13) 

This  equation  is  an  empirical  fit  to  tabulated  values  presented  in  reference  22. 
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TABLE  I.-  TEST  CONDITIONS,  DATA,  AND  CALCULATIONS  FOR  MOLYBDENUM 
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Specimens  are  numbered  for  convenience  in  identification  and  these  numbers  do  not  indicate  the  order  in  which  specimens  were  tested, 
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TABLE  II.-  EXPERIMENTAL  DATA  AND  CALCULATED  ACTIVATION 
ENERGY  VALUES  FOR  OXIDATION  OF  THIN  FLAT  PLATES  OF 
MOLYBDENUM  IN  TUBE  FURNACE 


Wall 

temperature, 

Tw.  °K 

(a) 

Flow 

velocity, 

ue, 

m-sec"l 

(a) 

Average 

oxidation 

rate,  kav, 
kg-m"2_sec"  1 

(a) 

Specimen 

dimensions, 

cm 

(a) 

Activation 

energy, 

Ea> 

J-(kg-mole)-'*' 

(b) 

Sei 

'ies  1;  air;  volumetric  flow  rate,  0.57  x  10"^  m3-sec“l 

1322 

1.76 

2.96  x  10'3 

2.49  x  1.14  x  0.25 

8.326  x  107 

1369 

1.83 

3.10 

2.49  X  1.14  X  0.25 

8.496 

1425 

1.90 

3.32 

2.49  X  1.14  X  0.25 

8.643 

1486 

1.98 

4.05 

2.49  X  1.14  X  0.25 

8.296 

1514 

2.01 

4.23 

2.49  X  1.14  X  0.25 

8.275 

Series  2;  air;  volumetric  flow  rate,  1.09  x  10"4  m^-sec--*- 

1314 

3.35 

3.25  X  10'3 

2.46  X  1.14  X  0.25 

8.459  x  107 

1347 

3.44 

3.18 

2.46  x  1.14  x  0.25 

8.717 

1386 

3.54 

3.76 

2.46  x  1.14  x  0.25 

8.589 

1436 

3.67 

4.40 

2.46  x  1.14  x  0.25 

8.489 

1450 

3.71 

4.86 

2.46  X  1.14  x  0.25 

8.268 

1478 

3.78 

4.91 

2.46  X  1.14  x  0.25 

8.396 

aFrom  reference  2. 


^Calculated  by  equation  (16). 
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TABLE  II.-  EXPERIMENTAL  DATA  AND  CALCULATED  ACTIVATION 


ENERGY  VALUES  FOR  OXIDATION  OF  THIN  FLAT  PLATES  OF 
MOLYBDENUM  IN  TUBE  FURNACE  -  Concluded 


Wall 

temperature, 
Tw,  °K 

(a) 

Flow 

velocity, 

ue> 

m-sec-l 

(a) 

Average 

oxidation 

rate,  kav, 

kg-m-2-sec"  1 

(a) 

Specimen 

dimensions, 

cm 

(a) 

Activation 

energy, 

Ea> 

J- (kg- mole)”  * 
(b) 

Series  3;  air;  volumetric  flow  rate,  1.49  x  10"4  m^-sec--*- 

1300 

4.55 

3.20  x  IO-3 

2.46  X  1.14  X  0.25 

8.526  X  107 

1350 

4.72 

3.57 

2.46  x  1.14  X  0.25 

8.640 

1369 

4.79 

3.52 

2.46  X  1.14  X  0.25 

8.787 

1400 

4.90 

4.15 

2.46  X  1.17  X  0.25 

8.626 

1425 

4.99 

4.74 

2.46  X  1.14  X  0.25 

8.457 

1455 

5.10 

5.08 

2.46  X  1.14  X  0.25 

8.450 

Series  4; 

21.5%  O2,  78.5%  He;  volumetric  flow  rate,  1.34  X  10" 

■4  rn3-sec-^ 

1306 

4.12 

4.42  x  IO-3 

2.36  X  1.19  X  0.25 

8.373  x  107 

1329 

4.19 

4.71 

2.46  X  1.35  X  0.25 

8.401 

1362 

4.29 

5.18 

2.39  X  1.14  X  0.25 

8.440 

1396 

4.40 

5.84 

2.39  X  1.14  x  0.25 

8.415 

1478 

4.66 

7.13 

2.41  X  1.14  X  0.25 

8.452 

Series  5; 

13.6%  O2,  86.4%  He;  volumetric  flow  rate,  2.07  X  10' 

•4  m3_sec-l 

1433 

6.97 

3.47  X  10“3 

2.41  X  1.14  X  0.25 

8.873  X  107 

1462 

7.11 

3.81 

2.39  X  1.14  x  0.25 

8.877 

1539 

7.49 

5.35 

2.39  X  1.14  X  0.25 

8.594 

1562 

7.60 

5.18 

2.46  X  1.17  X  0.25 

8.780 

1610 

7,83 

6.59 

2.46  X  1.17  x  0.25 

8.385 

1623 

7.90 

7.06 

2.41  X  1.14  X  0.25 

8.259 

1644 

8.00 

7.38 

2.41  X  1.14  X  0.25 

8.215 

aFrom  reference  2. 
^Calculated  by  equation  (16). 
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TABLE  III.-  COMPARISON  OF  CALCULATED  VALUES  OF  Ea 


Series 

Gas  composition 

Ea  from 
ref.  2, 

J- (kg- mole)- * 

Ea  from 
eq.  (16), 
J-(kg-mole)-'*' 

1 

Air 

8.53  X  107 

8.42  X  107 

2 

Air 

9.00 

8.49 

3 

Air 

9.16 

8.58 

4 

21.5%  02,  78.5%  He 

8.51 

8.42 

5 

13.6%  02,  86.4%  He 

9.37 

8.58 

Mean  value  of  Ea,  J-(kg-mole)-^  .  . 

8.91  X  107 

8.51  X  107 

Percent  standard  deviation . 

5.7 

2.1 

NASA- Langley,  1966  L-4245 
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"The  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  co?2tribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof 

— National  Aeronautics  and  Space  Act  of  1958 
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